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Abstract. An emissivity formula is derived using the generalised Fermi-Weizacker- 
Williams method of virtual photons which accounts for the recoil the charged particle 
experiences as it emits radiation. It is found that through this derivation the 
formula obtained by Sokolov et al using QED perturbation theory is recovered. The 
corrected emissivity formula is applied to nonlinear Thomson scattering scenarios in the 
transition from the classical to the quantum regime, for small values of the nonlinear 
quantum parameter y. Good agreement is found between this method and a QED 
probabilistic approach for scenarios where both are valid. In addition, signatures of 
the quantum corrections are identified and explored. 
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1. Introduction 

Increasing available laser power and high-energy electron beams (reaching few GeV; 
for e.g. m ei m obtained with plasma-based accelerators allow us to reach new 
regimes in laser-beam scattering in all-optical configurations H 5, [6j. At sufficiently 
high laser intensities, radiation damping becomes non-negligible and needs to be taken 
into account in the electron dynamics through the radiation reaction force. For even 
greater intensities or for extremely relativistic electrons and moderate intensities, the 
energy loss in the emission of a single photon can be so big that the trajectory is no longer 
continuous and radiation reaction needs to be described with quantum electrodynamics. 
The transition to the quantum regime can be parameterised by: 
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where E crit = m 2 c 3 /(eh) is the Schwinger critical held, F IW is the electromagnetic 
wave field-strength tensor and jf is the momentum four-vector. This parameter is 
approximately equal to the ratio between the held amplitude in the rest frame of 
the electron and the Schwinger critical held for very relativistic electrons. When 
X < 1 emission occurs through Compton scattering and, as the emitted energy becomes 
comparable to the electron energy, the trajectory gets stochastic; at x ~ 1 processes 
such as pair production become signihcant. However, for x “C 1, the trajectory of the 
electron can still be described classically if one accounts for the radiation reaction force. 
The emission process is Thomson scattering, the classical analog of Compton scattering. 
In this paper, the dehnition for nonlinear Thomson/Compton scattering from ref. [7J is 
used. 

Another important parameter in describing the nature of electron-laser scattering 
is the normalised vector potential of the laser oq. For very intense lasers, a () W 1 the 
electron dynamics in the laser is nonlinear and radiation is emitted at multiple harmonics 
(nonlinear Thomson [ 8 ; i9]/Compton H). Tl| scattering). In this work, the regime of 
radiation emission with x -C 1 to < 1 and ao 1 and 7 1 will be studied, where 7 is 

the Lorentz factor of the electron. In it, radiation damping is non-negligible but can be 
described classically and the higher harmonics of nonlinear Thomson scattering extend 
to energies approaching the energy of the electron, such that the classical emissivity 
formula derived from the Lienard-Wiechert potentials may not be valid anymore. The 
question then arises of how to model radiation emission in this scenario. 

Several techniques have been used, such as semiclassical calculations danann, 
QED perturbation theory m the introduction of functions that correct the equation of 
motion and the radiated power/spectrum [18, HI 17J and Monte-Carlo methods, either 
based on the cross sections for Compton scattering puns] or on the emission probability 
function and spectrum of synchrotron radiation combined [201 EH EH ESI EH]. 

In this paper we show how to derive a quantum corrected emissivity formula 
for arbitrary observation directions using the generalised method of Fermi-Weizacker- 
Williams, first developed by Lieu & Axford [2Bj. This formula is implemented in the 
post-processing radiation diagnostic code JRad EU. It is then used to explore the 
changes that occur in the Thomson scattering spectrum in the transition from the 
classical to the quantum regime. 

This paper is structured as follows: in section I, a quantum corrected emissivity 
formula is derived. In section II, comparisons are shown between the spectrum computed 
with this technique and the spectra obtained through a QED probabilistic approach 
for synchrotron radiation. In section III, we explore nonlinear Thomson scattering in the 
transition from the classical to the quantum regime with JRad and trajectories obtained 
from the integration of the equation of motion with the radiation reaction force [35J. 
We show that under certain conditions a unique signature of the quantum corrections 
is observed, which is not seen if only the radiation damping is taken into account in the 
trajectory. In section IV we state the main conclusions of the paper. 
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2. Radiation emissivity with quantum corrections 


The method of Lieu & Axford consists of an extension of the method of virtual 
photons (FWW method) [26]. This method was originally applied to problems such as 
bremsstrahlung radiation, where the electron was considered to be in uniform motion. 
The application of the FWW method to scenarios such as an electron gyrating in a 
uniform and static magnetic field is not possible, since there is no single rest frame 
for the electron. Lieu & Axford have proposed a method to overcome this issue [26j. 
Their solution relies on splitting the electron trajectory into a series of infinitesimal 
segments in which the electron velocity is approximately constant and therefore the 
method of virtual photons can be applied. They then derive the classical results of 
synchrotron radiation by determining the Thomson scattering spectrum in this series of 
instantaneous rest frames, transforming the emissivities back to the laboratory and then 
adding them coherently. Quantum corrections associated with the recoil of the electron 
during the emission process are added by replacing the Thomson scattering cross section 
by the Compton cross section and oj by ui/r], where // is the Compton shift. Though 
Lieu & Axford have extended their work to three-dimensional scenarios [2_8j, their 3D 
emissivity formula only reduces (in the absence of quantum corrections) to the classical 
result for certain angles of observation. 

In the following, it is shown that by introducing an additional generalisation to the 
method of Lieu & Axford, it is possible to obtain a quantum corrected emissivity formula 
which reduces to the three-dimensional classical emissivity in the limit of negligible 
Compton shift for arbitrary directions of observation. For the sake of completeness, 
the essential steps in the derivation are presented here; for further details, the reader is 
invited to check ref. [261127]. 

Firstly, a coordinate transformation is performed in the laboratory frame such that 
the segment to be analyzed has its beginning at the origin. A Lorentz transformation 
is then performed to the instantaneous rest frame of the electron for this segment. In 
this frame, the component of the Poynting flux for the frequency cu', Sf oj') is given by: 



( 2 ) 


where c is the velocity of light, E' is the incident electric field, t' is time and the primed 


variables refer to quantities in the instantaneous rest frame. 

In this infinitesimally small segment the electron will remain non-relativistic and 
the electric field can be related to the electron velocity from the equation of motion 
mir'_ l = eE', where m and e are the mass and charge of the electron, and v' ± is the 
acceleration of the electron. Assuming that the field can be decomposed into a series of 
plane waves, the Poynting flux for the segment then becomes after integration in time: 



( 3 ) 


where At' is the time the electron takes to cross the infinitesimal segment. 
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The emissivity at a frequency u / (energy radiated at a frequency u / per unit of 
frequency and per unit of solid angle) in the infinitesimal segment of trajectory, A a', 
is given by the product of the incident radiation flux and the Compton cross section 
(■ da'/dft')c = [e 2 /{me 2 )} 2 rj' |e(, ut • e' n | 2 , where e' in and e' out are the polarisation vectors 
of the incident and outgoing photons. To account for the Compton shift corrections, u>' 
is replaced by to a//?/ and the associated substitution: 


d 2 I 


d 2 I 
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is made. The emissivity in the infinitesimal rest frame with the quantum correction is 
then given by: 
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To transform this result back to the laboratory frame it is useful to note that the 
dot product e' out ■ e' in , v' ± (t')At' and (for very relativistic particles) rf are invariants, 
as pointed out by Lieu & Axford [26]. To facilitate the calculations, these two 
quantities are determined in the rest frame of the particle but expressed in terms of 
the variables in the lab frame [26.J. The transformation of the time and frequency yields 
ioj't' = iu(t — n • vt/c) and when transforming back from the simplified reference frame 
to the original laboratory frame the term n • vt in the exponent yields n • r. 

To obtain the emissivity in the laboratory, the contributions to the total radiation 
spectrum from each segment A a need to be Lorentz transformed to this frame and 
added coherently, by summing them before taking the square of the modulus |26j in 
equation (|5j): 
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To facilitate the calculations, the direction of the velocity can be taken to be in the 
positive x direction and spherical coordinates will be used. A natural choice of unity 
vector in the observation direction for the outgoing radiation is then n' = ej./|e'|, i.e. 
n' = (cos 9', sin 6' sin 0', sin 9' cos ft), where O' is the angle between the wave vector of 
the outgoing radiation and the velocity and ft is the angle between the two transverse 
components of the observation direction, z and y. The polarisation vectors can then be 
given by the other unity vectors in the spherical coordinates triad, e' e and ei, which in 
terms of the laboratory frame angles gives: 
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The displacement in the rest frame can be expressed in terms of the laboratory 
quantities: 


(Ax', Ay', Az') = (q(n x - v)At, v y At, v z At), 
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where 7 is the Lorentz factor of the electron, v x , v y and v z are the components of the 
velocity of the electron in the laboratory, in the trajectory infinitesimal segment. 

Substituting the rest frame displacement and the outgoing radiation polarisation 


vectors into equation (| 5 [) the specific emissivities in the polarisation directions e' 


out 


and el..* , denoted by a, and o : 2 respectively, can be computed. The two contributions 


~out 


can be summed inside the modulus since they are perpendicular to each other and the 
complete emissivity formula is obtained: 
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This formula (which recovers the result of Sokolov et al [17]) reduces to the classical 
emissivity in the limit 77 —>■ 1, whereas the three-dimensional result of Lieu & Axford 
[28] did not when <p 0 . The emissivity obtained here is then applicable to an 
arbitrary observation direction. The changes in the calculations above which allowed 
this generalisation were two: setting the unit vector of observation direction n' to 
n' = e^/le'j. and extending the perpendicular displacement vector Ar' = v'±At' to 
three dimensions. We also note that in doing the integration by parts in the calculation 
of cui, the term 
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was neglected since the modulus of a complex exponential is < 1 and it is assumed that 
co 1 . 

Equation © can be implemented in post-processing radiation diagnostic codes 
such as JRad [31], which was the one used in this work. The particle trajectories can be 
obtained from particle-in-cell codes or from the integration of the equation of motion, 
including the radiation reaction force. 


3. Comparison with QED results 

As a test to the quantum corrected diagnostic (JRad-QC), the synchrotron radiation of 
an ultra-relativistic electron in an ultra-intense magnetic field was computed with both 
JRad-QC and with OSIRIS-QED j2H 23J- The latter calculates the QED probability of 
radiation emission and determines the photon energy to be emitted using a Monte-Carlo 
method and the QED synchrotron radiation spectrum. 

In this benchmark the spectrum of an electron with 7 = 200 subject to a magnetic 
field of B — 5.7 x 10 10 G (which corresponds to 7 = 0.26) is compared with the 
result from OSIRIS-QED without damping. The case without damping is presented 
for reference since the theoretical radiation spectrum is known, and therefore can easily 
be used for benchmarking purposes. In this case, the damping is turned off in OSIRIS- 
QED, hence the spectrum comparison tests solely whether the quantum correction in 
the model of Lieu & Axford reproduces the results obtained for photon production 
according to a QED calculation of nonlinear Compton scattering. 
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Figure 1. (left) Spectrum determined over a line perpendicular to the synchrotron 
trajectory plane determined from the JRad diagnostic with quantum corrections 
(JRad-QC). (right) The red line represents the synchrotron spectrum obtained with 
JRad-QC and the black line is the result obtained with OSIRIS-QED. 


For this benchmark, an analytical trajectory was produced according to the classical 
equations of motion for an electron in a static magnetic field (see for e.g. [31)j). The 
post-processing diagnostic with quantum corrections was used to obtain the emitted 
spectrum as a function of solid angle along a line perpendicular to the trajectory plane, 
at a distance of about 150 c/uy from the gyration circumference (the precise value is not 
relevant as long as it is in the far field), where uy = 1.88 x 10 15 rad/s is the normalisation 
frequency, equivalent to E = hoj r ~ 1.24 eV. This spectra was then integrated in solid 
angle to obtain dE/doj (seen in red in Figure [l|. The result shows excellent agreement 
between the OSIRIS-QED and the JRAD-FWW spectra (which had been shown to be in 
agreement with the theoretical result for synchrotron in the QED regime). We observe 
that comparisons with QED calculations need to be done with careful consideration 
since the method is only valid for ultra relativistic particles and sufficiently high fields 

Isa mi ng. 

4. Transition from the classical to the quantum regime 

To explore the nonlinear Thomson scattering spectrum in the transition to the quantum 
regime, trajectories from ultra-relativistic electrons colliding with plane waves and laser 
pulses were produced by numerical integration of the equation of motion with a Runge- 
Kutta of 4 th order [35 1. The radiation diagnostic JRad-QC / JRad was then used to 
compute the radiation spectra at line detectors in the far-field region with/without 
quantum corrections. 

To investigate when quantum effects start to affect significantly the radiation 
spectrum a series of spectra has been computed for the case of a circularly polarised plane 
wave with normalised vector potential of a 0 = 0.88, 1.77, 7.07 scattering off an electron 
with 7 = 10 3 (travelling in the +x\ direction) in counter-propagating geometry. The 
plane wave is preceded by a rising ramp but the electron energy is ajusted so that 7 cx 10 3 
at the beginning of the flat part of the wave, which lasts for tfi at ~ 38 T 0 — 127 fs, where 









Modelling radiation emission in the transition from the classical to the quantum regime7 

T 0 is the wave period. The virtual detector is a line at £3 = 0 and the normalisation 
frequency is the same as the frequency of the monochromatic wave ujq = 1.88 x 10 15 
rad/s (which corresponds to 1 pm wavelength or E = ftwo ~ 1.24eV). 

Two effects are readily observed in the spectra in Figure |3j Firstly, in the case 
of a 0 = 0.88 there is a shift in the frequency of the harmonics. Secondly, in all the 
cases, a reduction of the radiated energy at higher frequencies is observed, as expected 
[14, 43]. The energy lost by the electron over the interaction time for the case of 
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Figure 2. Spectra of the scattering of a circularly polarised plane wave with ao = 0.88 
(1), a 0 = 1.77 (2) and a 0 = 7.07 (3). Plots labeled with (b) refer to the calculation with 
quantum correction (JRad-QC) and those labeled with (a) to the classical emissivity 
calculation (JRad). The white lines represent the spectra integrated in solid angle 
assuming cylindrical symmetry. 

a 0 = 7.07 was about 21%, which gives an average energy loss rate of 0.55% per period of 
oscillation. This slow energy loss rate is the reason behind the almost symmetric profile 
in the spectra, which would not be expected if the electron would change its energy 
by a significant amount during an oscillation period. It is then possible to estimate 
the total energy captured by the detector by doing an integration in solid angle and 
in photon energy, assuming cylindrical symmetry around the direction defined by the 
initial momentum direction of the electron, in this example the line in the +aq direction, 
at £2 = 0 and .1:3 = 0. Doing so, one obtains 271 keV from the integration of the detector 
without quantum corrections, 104 keV for the detector with quantum corrections, and 
109 keV from the direct measure of the energy loss by the particle in its trajectory. 
The energy captured in the detector with the quantum corrections is very close to the 
value measured in the track, which further supports the validity of the formula for the 
quantum corrected emissivity. 
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Regarding the other values of ao investigated, the energy captured in the detector 
with quantum corrections is also closer to the energy difference measured in the track 
but the difference is smaller, which could be anticipated from the fact that the observed 
differences in the spectrum are also less significant than in the highest ao case. 


5. Laser pulse scattering 

The scattering of a linearly polarised laser pulse from an electron with 7 = 1030 
(travelling in the +x\ direction) in counter-propagating geometry was also explored, in a 
configuration similar to ref. [lj. In the scenario studied here, the laser pulse has a 1 /j m 
wavelength, 26.5 fs duration and peak normalised vector potential of a 0 = 10, 20, 30 
and is linearly polarised in the x 2 direction. The spectra computed for a line positioned 
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Figure 3. Spectra of the scattering of a linearly polarised laser pulse with ao = 10 
(1), ao = 20 (2) and ao = 30 (3). Plots labeled with (b) refer to the calculation with 
quantum correction (JRad-QC) and those labeled with (a) to the classical emissivity 
calculation (JRad). 

at X\ = 10 1 c/ujo and x 3 = 0 e/cq, with and without quantum corrections are depicted 
in Figure [3| From the results, it can be seen that as the peak ao increases the shape of 
the spectrum with quantum corrections changes more and more significantly, specially 
at higher angles of observation (higher x 2 values in the detector axis). In addition, at 
the highest ao, spike structures are clearly seen. These are correlated with the maxima 
in the transverse momenta, where the derivative of 7 is zero (see Figure [4]), and exist 
also in the classical case. However, the quantum corrections lead to a stronger decrease 
in radiated energy between the spikes, making them appear more clearly. 

The spectrum changes are more complex and seem more dependent in angle in 
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Figure 4. Evolution of the electron transverse momentum (p2) (left top pane) while 
interacting with a laser pulse with ao = 30 and of the position in the detector to 
which it points to (bottom left pane). Line spectrum from the scattering of the laser 
pulse off the electron (right pane). The dashed lines show the peaks in the spectrum 
originate from points in the trajectory where d^/dt = 0, which correspond to peaks in 
the transverse momentum. 


the laser pulse scattering setup compared with the plane wave scenario. This will be 
explored in more detail in a future publication. 

6. Conclusions 

In this paper it was shown that by extending the generalised FWW method of Lieu 
& Axford, the emissivity with quantum corrections due to the electron recoil can be 
derived, which had been obtained by a different approach by Sokolov et al HZ], through 
QED perturbation calculations. 

The quantum corrected emissivity was implemented in the numerical diagnostic 
JRad EU and used to determine the nonlinear scattering spectrum of relativistic 
electrons with ultra-intense laser pulses and plane waves. It was found that only when 
the quantum corrections are introduced does the energy captured in the virtual detector 
spectrum become consistent with the energy loss by the particle (as measured from the 
integration of its equation of motion). 

The analysis of the scattering of ultra-intense linearly polarised laser pulses from 
relativistic electrons revealed that the changes in the transverse momenta and energy of 
the electron during the interaction led to a complex spectrum shape. For larger angles 
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of observation, spike features are observed in the spatially resolved radiation spectrum, 
surrounded by regions of lower radiation emission. The shape changes reflect the changes 
in the Compton shift during the interaction, which were much more significant in the 
laser pulse scenario than in the plane wave case for the parameters used. This can be 
attributed to a combination of stronger radiation damping and bigger changes in the 
field amplitude experienced by the electron in the laser pulse scattering scenario with 
higher peak a 0 as compared to the plane wave cases. 
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